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Abstract 

Using the Fokker-Planck equation describing the evolution of the transmission 
eigenvalues for Dyson's Brownian motion ensemble, we calculate the magne- 
toconductance of a ballistic chaotic dot in in the crossover regime from the 
orthogonal to the unitary symmetry. The correlation functions of the trans- 
mission eigenvalues are expressed in terms of quaternion determinants for 
arbitrary number N of scattering channels. The corresponding average, vari- 
ance and autocorrelation function of the magnetoconductance are given as 
a function of the Brownian motion time t. A microscopic derivation of this 
S-Brownian motion approach is discussed and t is related to the applied flux. 
This exactly solvable random matrix model yields the right expression for the 
suppression of the weak localization corrections in the large A^-limit and for 
small applied fluxes. An appropriate rescaling of t could extend its validity to 
larger magnetic fluxes for the averages, but not for the correlation functions. 
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I. INTRODUCTION 



Ballistic transport through quantum dots of various shapes are the subject of many 
theoretical and experimental investigations. One of the fundamental issues is to determine 
in a conductance measurement the signature of "quantum chaos", i. e. of the quantum 
analog of a classically chaotic dynamics. Usually, such dots are made with two dimensional 
semiconductor nanostructures (quantum billiards) connected to two electron reservoirs, and 
one measures the magnetic field dependence of the conductance at very low temperature 
(quantum coherent transport). The question is to describe the statistical properties of 
the recorded magneto-fingerprints for different shapes of the dot. We restrict ourselves to 
fully chaotic systems where we hope to state universal results, leaving aside the integrable 
systems. The usual theoretical approaches concentrate on microscopic models, using 
numerical calculations and semiclassical approximations, or more or less straightforward 
random matrix models, possibly combined with supersymmetric methods. The semiclassical 
approach requires in the final evaluation the so-called diagonal approximation, which is 
problematic |P for the total weak localization correction. Supersymmetry and random 
matrix model have also been recently used |^ to describe the ensemble averaged behavior 
of the magneto conductance. In Ref. 0] the Hamiltonian of the dot was modeled by a 
Pandey-Mehta random matrix Hamiltonian 0, but the used method leads to substantial 
technical problems and the averaged magneto conductance is expressed in terms of a rather 
complicated integral expression which can be evaluated by numerical means. 
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Another random matrix approach for the quantum biUiards was used in Refs. |TO|jTl[ 
models directly the scattering matrix S of the dot instead of its Hamiltonian. There are 
reasons |l2| , p!3| to assume that S is suitably described by one of Dyson's classical circular 
ensembles (COE-CUE-CSE). In Refs. |1I0| . |TT[| , the associated transport properties have 
been calculated: e. g. for the circular orthogonal ensemble (COE) and for circular unitary 
ensemble (CUE) corresponding to chaotic dots without or with a large enough applied 
magnetic field respectively. We propose a generalization of this approach in terms of a 
continuous crossover from COE to CUE, that is described by one of Dyson's Brownian 
motion ensembles |]TB|. The idea is to consider a Brownian motion with a fictitious time t 
that leads to a diffusion in ^-matrix space. In the limit t ^ oo the probability distribution 
diffuses to a stationary solution that corresponds in our application to the circular unitary 
ensemble. Choosing as initial condition for t = the circular orthogonal ensemble, we have 
a model for the magnetoconductance of a chaotic dot, at least as a function of the Brownian 
motion time t. Before completion of this manuscript, we received a preprint by Jochen 
Rau [|16| where a very qualitative analysis of this model is presented with essentially similar 
conclusions than our detailed study. One can also mention that this crossover ensemble 
was already considered in Refs. [lT7| , p!8| where the correlation functions for the scattering 
phase shifts 6i were calculated. The decisive difference here is that we are interested in the 
transport properties, i. e. in the statistics of the transmission eigenvalues Tj of tt\ when S 
is parametrized by 



The Brownian motion leads to a Fokker-Planck equation for the probability distribution of 



the variables Tj which have been first derived in Ref. [1^. This Fokker-Planck equation 
can be mapped onto a Schrodinger equation of a one dimensional system of free Fermions 
and solved for an arbitrary initial condition. The corresponding propagator was already 



calculated in Ref. where the crossover of two decoupled CUEs of dimension to one 
CUE of dimension 2N was considered. Here we use another initial condition that is just 
the joint probability distribution for the variables Tj of the COE-case. The time dependent 
solution of the Fokker-Planck equation is then the joint probability distribution for the 
crossover ensemble COE CUE. 

This paper is subdivided in two main parts (Sections y and |T|). Section ^ is completely 
devoted to an exact solution of the Brownian motion model as it stands. We give any 
/c-point correlation function for the Tj at one particular time t. In addition, we give also 
the most simple correlation functions between two different time values. The technical 
method to obtain these results is based on a formulation in terms of quaternion determinants 
20| and skew orthogonal polynomials. The final results for the correlation functions are 



given in terms of Legendre Polynomials and can be easily evaluated by numerical means. 
The averages for the conductance and the autocorrelation of its variance (between different 
Brownian motion times) are exactly calculated for an arbitrary value of the channel number. 
All results of Section |T| are exact consequences of the S'-Brownian motion ensemble applied 
to the COE — * CUE crossover. 

It remains to express the Brownian motion time t in terms of a real physical parameter 



such as the applied magnetic flux $. This issue is considered in detail in Section III. We 
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first discuss the relation of the Pandey-Mehta Hamiltonian 0] with finite disordered systems 
(in the diffusive regime) or with quantum bilhards (in the chaotic balhstic regime). Both 



cases have already been considered in Refs. pT|-p3| and we give a simple treatment that 



accounts for both cases and generalizes the known results to the quantum billiard case. In 
subsection [11 B| we briefiy recall the different magnetic field scales in a quantum billiard 
and the effects of the Landau quantization for very strong fields. In subsection [III Q , a 
microscopic justification^ of the Brownian motion ensemble is sketched, assuming a crossover 
Pandey-Mehta Hamiltonian for H and using a well known expression of 5* in terms of H. 
We obtain an explicit expression that relates the Brownian motion time with the Pandey- 
Mehta parameter a and therefore with the magnetic fiux. Using this translation we find a 
Gaussian instead of a Lorentzian behavior (as given in Ref. [§) for the suppression of the 
weak localization correction to the averaged conductance, when a fiux $ is applied. One 
finds a certain characteristic fiux $^7 separating two regimes. For a small fiux $ < $^7 our 
result matches precisely that of Ref. |^ . For a large fiux $ 3> $c the discrepancy between 
the Gaussian and Lorentzian becomes more important. This indicates that the Brownian 
motion model is not precisely equivalent to the model used in [Q for large times of the 
Brownian motion. The difficulty to obtain directly from our model the lorentzian shape 
for the weak localization suppression and the lorentzian square behavior for the correlation 
functions is discussed in conclusion. 



II. 5-BROWNIAN MOTION ENSEMBLE AND COE TO CUE CROSSOVER 

In Ref. |TP| we used the original Dyson Brownian motion ensemble for the scattering 



matrix to obtain a Fokker-Planck equation for the transmission eigenvalues Tj. This is 
suitable for having the conductance g oi & scatterer connected by A^-channel contacts to 
the electron reservoirs, since g = J2iLi Ti. The stationary solution of the Fokker-Planck 
equation is given by the transmission eigenvalue distribution of the circular unitary ensemble 
describing a chaotic dot without time reversal invariance. The chosen initial condition 
for the Brownian motion is the COE-distribution for the Tj. Concerning the definition of 
the Brownian motion, the derivation of the Fokker-Planck equation for the transmission 
eigenvalues and its solution for an arbitrary initial condition, we refer the reader to Ref. 
I9| . It is convenient to use the same coordinates Xi = 1 — 2Ti of [0 , where the Tj are the 



eigenvalues of tt^. 



A. Joint probability distribution of the transmission eigenvalues 

The joint probability distribution p{x, t) for the Brownian motion ensemble can be ex- 
pressed as follows 



^ We mention that in Ref. |18| a semiclassical and numerical justification for the validity of the 
Brownian motion ensemble has been given as far as the scattering phase shifts correlations are 
concerned. 
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p{x,t) = J d^x p{x)G{x,x;t) , (2.1) 
where the many particle Green's function is given by 

Gix, X- 1) = -^P{x)-' p{x) det(^(x„ X,- 1)) e-^^* , (2.2) 

N 

P{^) = X{{^^-^3) , (2-3) 

i>j 
oo 

g{xi, xf, t) = E o (1 + 2^) Pn{xi) Pn{ij) e-""* . (2.4) 

n=0 ^ 

Pn{xi) is the Legendre polynomial of degree n, £„=(! + 2?t,)^ is the eigenvalue of the "one 
particle Hamiltonian" [O] and the constant Cn is given by the sum = —J2n=o^n- The 



n=0 ^n- 

function p{x) = p{x, 0) is the initial condition. In our case it is given by the joint probability 
distribution for the orthogonal circular ensemble which in terms of the variables xj reads as 

SHI 



N N 

p{x) OC Jl = n ~ ^jl • (2-5) 



i=i y'2(l - Xi) i>j 

The integral in Eq. (|2.1| ) can be evaluated by the method of integration over alternate 
variables ||2^. As usually in such situations, we consider now the case of an even number 
of channels. The complications that arise from the case of odd N will be discussed in 
Appendix |B[ Our final results of this Section are valid for all values of A^. The result for 
the joint probability distribution can be expressed in terms of the so-called Pfaffian 

N 

p{x, t) OC n l^i - I \/ det {H{xi,Xj;t) e"*^^* (2.6) 

i>j 

where H{xi, xj; t) is an antisymmetric function given by the double integral 

H{xi,X2;t) = ]- f dyi [ dy2 . ^^^^ 9{xi,yi;t) g{x2,y2;t) (2.7) 

2 J -I J -I ^(l-y2)(l-yi) 

where e{u) = +1, 0, —1 according to the cases u > 0, =0, < 0. In the next subsection 
we will present a formulation in terms of quaternion determinants ||25[ that enables us 
to evaluate all types of correlation functions and to get exact expressions for the average 
conductance and the conductance fluctuations. 

We first briefly discuss the implications of Eq. (|2.6|) for the "level" repulsion in a more 
qualitative way. For this discussion the angles ipj with Xj = sin^ ipj are more appropriate. 
The one particle Green's function g{xi,Xj;t) fulfills in terms of the (fj a time dependent 
Schrodinger equation [|T^ with a particular potential. In the (extreme) short time limit one 
can neglect the effect of this potential and the Green's function can be approximated by the 
free "diffusion propagator", i.e. 

g{siia^{^i), sin^i^j); t) ~ , ^ ^ ^ exp - fj^) ■ (2-8) 
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Using this expression we can evaluate the integral in Eq. ( |2.7| ) and we obtain for the joint 
probability distribution for the variables (pi in the limit t — the expression 



~Fi(¥;)^det (erf -</.,)/v/4t)) (2.9) 

where Fi{ip) is the joint probability distribution for the orthogonal case, i.e. for exactly 
t = 0. Pandey and Mehta ^ obtained for the eigenvalue distribution of the random matrix 
Hamiltonian described by Eq. ( |3.1| ) a very similar expression if the (pj are identified with 
the energy levels of and if one relates the symmetry breaking parameters a and t in a 
suitable way. The expression ( |2.9| ) shows that the crossover from COE to CUE has an effect 
for arbitrarily small values of t > if the difference \ipi — (pj\ is sufficiently small. Let us 
first consider the situation where all ipi are well separated, i.e. \'Pi — 'Pj\ ^ -\/4t. In this case 
the error function in Eq. ( p.9|) can be replaced by sign((y9j — ipj) and the determinant is just 
one. The probability distribution then equals the distribution Fi{ip) for the COE. On the 



other hand the determinant that appears in (2^) is an even function of — (fj and vanishes 



if one of these differences is zero. If for example two values of the (pi are very close, e.g. 
|v?i — V'2| ^ "\/4t, the square root of the determinant becomes proportional to the difference 
— P'2\/V4i which leads together with the factor Fi((p) to the quadratic level repulsion of 
the unitary case. This behavior means that the COE to CUE crossover for the transmission 
eigenvalue correlations appears at arbitrarily small values for the time parameter t on scales 
smaller than -y/it. The same effect is also observed for the energy levels Ea [pi and for the 



scattering phase shifts 9i [17 



B. Quaternion formulation 

In this subsection, we use the technic of quaternions, that was introduced in the theory 
of random matrices in Ref. |2^, to calculate the m-point correlation functions given by the 
usual definition 

Rm{Xi,. . . ,Xm;t) = — -/ dXm+1--- dXNp{x,t) . (2.10) 

[Jy — my. J-i J-i 
We define for two arbitrary functions /, g the antisymmetric scalar product 

<f,9 >^*^= J'^dy, J'^dy, H{yuy2;t) f{yi) giy^) (2.11) 

where H{yi, 7/2; t) is just given by Eq. ( p.7|) . Suppose that the number of channels is even 
and that we have a set of skew orthogonal polynominals q':^\x) of degree n = 0, . . . , — 1 
with respect to the scalar product ( p.ll| ), i.e. 



<gi^g^^>^*^=^n™ (2.12) 

where is an antisymmetric matrix with Z2n,2n+i = —Z2n+i,2n = 1 and all other Znm = 0. 
The superscript means that both the polynomials and the scalar product depend explicitly 
on the time parameter t. At the moment, we do not consider explicit expressions for the 
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q^\x) and the following method remains general. In addition, we introduce the functions 

QW(x) via 

Qli\y2) = j\dy,H{y,,y^-t) q^\y,) (2.13) 

which are dual to the q^\x) : 

l'^dyQii\y)q^^\y) = Z^m • (2.14) 
We define furthermore functions of two variables x, y by 

N-l 

Kpp{x,y;t)= Y: Fi'\x)Z„^F(^\y) (2.15) 



n,m=0 



where F and F stand for one of the symbols q or Q. Hence, Eq. ( |2.15|) defines four different 
functions Kqq, KqQ, Kqq and Kqq. From the definition of the matrix Znm and Eqs. ( p.l3 - 
2.15| ) we can immediately state the following properties 

Kpp{x,y;t) = -Kpp{y,x;t) , (2.16) 
- J ^dx KqQ{x,x;t) = J ^dx KQq{x,x;t) = N , (2.17) 

J^^ dy Kpg^x, y; t) KQp{y, z) = - J^^ dy Kfq{x, y; t) K^p{y, z) = Kpp{x, z; t) , (2.18) 



and 



I'^dy H{y,x-t) KqFiy,z-t) = KQp{x,z-t). (2.19) 



From now on, we strongly refer to the notations and results of Ref. |]25| concerning quaternion 
properties. We introduce the 2 x 2-matrix (or quaternion) function a{x, y; t) by 

<^i^^y^t) = ( m Kqqix,y;t) \ 

^ ^ \H{x,y;t)-KQQ{x,y;t)KQq{x,y;t) ^ ^ 



that fulfills the three properties 



1 



-1 



dx a{x,x;t) = I 1 (2.21) 



1 



dy a{x,y;t) a{y, z;t) = (j{x, z;t) + T (t{x, z;t) — a{x, z;t) T , (2.22) 
p{x, t) = const. ■ QDet {{a{xi, Xj)i<ij<N) (2.23) 



where t = |Q and QDet(- ■ ■) is the quaternion-determinant defined in Ref. p5|. The 
verification of Eqs. (|2.21|) and (|2.22| ) is a straightforward calculation applying the properties 
( p.l6| - |2TT9| ). The proof of Eq. ( |2.23| ) is not so obvious and can be obtained as follows. First, 
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we note that ( p.l6|) and the antisymmetry of H{x,y;t) imply that the quaternion-matrix 
a{xi,Xj;t) is self-dual, i.e. 

a{xi,Xj;t) = a{xj,Xi;t) (2.24) 

where 



a h \ _ ( d —h 
c d I I — c a 



(2.25) 



is the quaternion adjoint in the notation of Ref. p5|. As in Ref. I^Sf we denote with 
A((j(xj, Xj\ t)) the conventional 2N x 2A^-matrix that contains the quaternions as 2 x 2 blocks. 
From theorem 2 of ||2^ we get the identity (QDet((j(a;j, x^; t)))^ = (lQi{A{a{xi^Xj;t))) that 
relates the conventional with the quaternion-determinant for self-dual matrices. We consider 
therefore the matrix 

A(a(x- x 't))= \ ~KqQ{^i-:^j\'^) Kgg{xi,Xj;t) \ 

^' ^' \H{xi,Xj;t)-KQQ{xi,Xj;t) KQq{xi, Xj;t) J 



qix.f 0\ Z 0\ -Q{x,) qix, 
Q{x,f 1 [O 1 [H{xi,x,;t) 



(2.26) 



Here we have used a quasi block notation where 1 and Z stand for the N x N unit matrix 
or the matrix with entries Znm respectively. q{x) (or Q{x)) is a column vector with entries 
QnK^) i'^^QnK^))- 5'(a;)^ and (5(a;)"^ denote the corresponding row vectors. The determinant 
of ( ^.261 ) is immediately calculated 

det{A{a{xi,Xj;t))) = det (^g^*^(xi))^ det{H{xi, Xj;t)) ~ (p(x,t))^ (2.27) 

where we have expressed the determinant of the skew orthogonal polynomials as a Vander- 
mond determinant. Eq. ( |2.27| ) completes the proof of the property ( |2.23| ). 

We can now use ( |2.21[ - 12.23| ) and apply the theorem 4 of Ref. [|25|. Then, we find directly 
the m-point correlation functions as quaternion-determinants 

Rm{xi, . . . ,x,n;t) = ClDet{a{xi,Xj;t)i<ij<ra) . (2.28) 
In addition, the constant in ( p.23|) must have the value 1/N\. 



C. Conductance correlation functions and averages. 

A more explicit evaluation of Eq. ( p.28|) for the m-point correlation functions requires 
the knowledge of the skew orthogonal polynomials g^*^ (x) and of their dual functions Q^^ (x) 
given by ( |2.13| ). The calculation of Appendix ^ shows that the skew orthogonal polynomials 
are not unique and one has to specify a certain choice. Of course, the final results of the 
last subsection do not depend on the particular choice. A possible choice is 

gS(a;) = i(l + 4n) P^n^x) e^-* + 5^(2;) , 
q2l+M = 1(3 + 4n) P2n+i{x) e^-+^* + 4*)(x) (2.29) 
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where Pnix) are the Legendre polynomials and 



2n-l 



^(-l)™i(l + 2m)P„(x) e 



(2.30) 



m=0 



The complete calculation that yields Eq. ( 2.29 ) is done in Appendix]^. Eq. (|2.29 ) can be 
inverted with the result 



x] 



x] 



i(3 + An) P2n+i{x) e=-+^* = g^li(x) - ^{x) 
where now the sum s^^ (x) is expressed in terms of the g^*-* (x) 

2n-l 



sii\x) 



X 



(2.31) 



(2.32) 



m=0 



We consider now Eq. (|2.31| ) for t = and replace in the definition of g{x,y;t) (compare 
Eq. (p.4|) ) one of the Legendre polynomials Pn{y) with the expansion in terms of the qn\y)- 
Then the integral in Eq. ( p.7|) can be done because the q^Ky) are just the skew orthogonal 
polynomials of the scalar product (|A1|) . After some algebra with the sums (it is useful to 
consider these type of operations as matrix multiplications with suitable defined matrices) 
we arrive at the expansion of the kernel H{xi, X2] t) 



H{xi,X2]t) = ^ e 

k,n=0 



-(£n+Sk)t 



sign(A; - n) Pn{xi) Pk{x2) 



(2.33) 



Our first application of ( ^.33[ ) is the calculation of the dual functions Q^^{x) that appear in 
the Eq. (|2T3| ), 



Q 



where we have introduced the sum 



Q^l = P2n+i{x) e-^-+^* + 5W( 
^li = -P2n(x)e--* + ^W(x) 



X 



oo 

E ' 

m=2n+2 



PM 



(2.34) 



(2.35) 



We have now to insert the expansions (^^), (^M ) and (lOil ) in the Eqs. ( ^^201) , (^^M 
in order to get the m-point correlation functions. In the following, we state the result for 
the 1- and the 2- point functions. First, we emphasize that this quaternion description, as 
it is shown in the last subsection, is valid for an even value of but the formulas below 
are valid for the odd A^ case, too. In Appendix |B|, we discuss briefly how the complications 
for odd A^ can be treated. Second, we give also an expression for the correlation function 
Ri^i{x, t;y,t + t) that describes the probability density of finding x at time t and y at time 
t + T. Its formal definition is given in and we calculated this function using a standard 



method of functional derivatives described in the book of Mehta 
For convenience, we introduce the following three functions 
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N-l 

St''\x,y) = Y: i(l + 2n)P„(x)P„(i/)e^"(*^-*^) 

n=0 

N-l oo 

+ ^(_l)iv-i-fc i(l + 2fc)P,(x)e^^*^ J2PMe-'"'' , (2.36) 

A:=0 n=N 
N-l 

Dp'''\x,y) = ^ sign(n-fc)i(l + 2n)(l + 2A;)(-ir+'=P„(x)Pfc(i/)e^"*^+^'=*^ , (2.37) 

n,k=0 
oo 

l'}^''''\x,y) = - sign{n-k)Pnix)Pkiy)e~'"'^-"^'' . (2.38) 

n,k=N 



Then the quaternion function ( p.20| ) becomes 



I^^'*\x,y) S}J'*^(y,x) 



and we eventually get for the correlation functions 

Ri{x;t) = S^j^''\x,x) , (2.40) 

R2{x,y-t) = R,{x-t) R,{y-,t) - Sp'\x,y) S'}^''\y,x) + D^^''\x,y) , (2.41) 
Pi,i (x, t; t + r) = Pi (x; t) Pi (y ; t + r) - (x, y) S^^^''^ (y , x) 

+D^^'^^\x, y) /g'*-^^)(x, y) + g{x, y- r) S'}^-^^''\y, x) . (2.42) 

For r = the function ( |2.42| ) fulfills the identity 

Pi,i(x, t; y, t) = R2{x, y; t) + 6{x - y) Ri{x; t) (2.43) 

where the last term accounts for the selfcorrelation of the eigenvalue x. In the limit ti = 
)f:2 — ^ oo in (|2.36|) only the first sum survives and the product ■ 1^ vanishes. Then Eqs. 
( [^.4U| - P7¥I| ) equal the 1- and 2-point correlation functions of the CUE case that are directly 
obtained with the Legendre polynomials as orthogonal polynomials [jlT|. The first corrections 
to this behavior for t 3> 1 are proportional e"*^^^""^^-!)* = e~^^* = e^*/*'^ where tc = 1/{8N) 
is the critical time for the Brownian motion, already given in |]19|. Of particular interest 
is the density for one channel which is just for N = 1 the distribution of the conductance 
(7 = T = (1 — x)/2 in units of 2e^ /h. The result ( |2.40| ) implies 

oo 

p(T; t) = 2Pi(l - 2T; t)\^^^ = ^ P„(l - 2T) e-^<^^^^' . (2.44) 

n=0 

The limiting cases are p(T; 0) = 1/{2\/T) and lim(^ooP(7'; t) = 1 and correspond to the 
results of the orthogonal and unitary cases respectively given in Refs. |lTO| , p]T| , which have 
been found to agree with microscopic numerical results |11|. Fig. 1 illustrates the crossover 
between these limiting cases, where p{T]t) is shown for t/tc = 0, 0.5, 1, 2, cxo. 

In the limit oo, it is possible to get a more explicit expression for the density. It 

is useful to express the density in terms of the angle variable cp, i.e. to consider pi{ip]t) = 
2sin(2(y9) Ri{cos{2ip) , t) . The most important contributions in the second sum of Eq. (|2.36|) 
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come from terms with indices near N, i.e. N — k N and n — N N . We can therefore 
hnearize Sj. — Sn — —8N{n — k) and replace the Legendre Polynomials by their asymptotic 
expression for large n, k. After some smoothing over the strongly oscillating contributions 
we get the expression 



Pi(<^;^) 



^-(n + J: e-(^'+i)('/'^) cos(2(2/ + 1)^)^ =^(n+ ^^"^^'/'^^ ^^^^'^^ 



71 



1=0 



cosh(2t/tc) — cos(4y9) 



which simplifies in the limit t ^ to 

2N 



p(^;0) = ^ + ^(5+M-5+(f-v')) 



(2.45) 



(2.46) 



where ■ ■) is the half-sided delta function that is only defined for positive arguments. In 
terms of the transmission eigenvalues T the density Ri{T; t) = 2 — 2T; t) reads 



Ri{T;t) 
with the limit 



1 



TT 



T(l -T) 



N + 



smh{t/tc 



2T) 



2(sinh^(t/t2)+4T(l-T))^ 



(2.47) 



^i(T;0) 



N 1 1 

, , +j{6^{T)-6^{l-T)) 

^ JT(1 - T) 4 



f2.48) 



We see that the asymptotic density for the variable (p is constant for the leading order in N. 



Eqs. (|2.45| - pl48|) contain also the next order, responsible for the so-called weak-localization 
correction of g. Eq. ( p.48| ) confirms the result of Ref. 



Some finite A^-results for the density are given in figures 2 and 3 (density for the trans- 
mission eigenvalue T with = 5). Fig. 2a and 2b compare the CUE- and COE-limits with 



the asymptotic density (cp. (|2.47|) ). One can find strange that the COE-density at finite 
N apparently fits better the asymptotic density than the CUE-density. This is due to the 
finite A^-oscillations on the scale of the "level-spacing" which are stronger for the CUE than 
for the COE. In the asymptotic results ( 2.47|) and ( 2.48 ) these oscillations are smoothed 
out. The crossover density interpolates between the two limits. Since the relative change 
for iV = 5 is not so strong as for = 1, in Fig. 3 only the difference Ri{T, t) — Ri{T, oo) 
is shown for t/tc = 0, 0.5, 1, 2, oo on a larger scale. In Fig. 4 a representative range for 
the 2-point correlation function i?2(To,T;t) for A^ = 5, t/tc = 10~^, 10~^, 10~^, 1, oo and 
To = 0.7, 0.6 <T< 0.8 is illustrated. One can very well see the effect already discussed at 
the end of subsection P^I A| that the COE CUE crossover happens at first for very small 
differences |T — To| iit/tc is small. All the curves shown in Fig. 1-4 are obtained by a direct 
numerical evaluation of the expressions ( p.40[ - pI^iD and ( |2.36| - p738| ). 



^The delta function at T = (or ip = 0) was not found in [l^, since the density was expressed in 
terms of the variable A = l/T— Iso that T = corresponds to A ^ oo, but this contribution is 
needed here to conserve the normalization. 
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The average of the conductance g = J2j Tj = and the autocorrelation of the 

fluctuation 6g = g — (g) can be expressed in terms of integrals with the densities ( p.40 - p^ ). 
These integrals can be evaluated by the recursion relation of the Legendre polynomials and 
their orthogonality relation with the result 

(^g(^t)) = -- , ^ , e-«^* , (2.49) 
2 2(2iV+l) ' ^ ' 

(5a(t) Sa(t + t)) = e^®^" I — + + ^^s(2N+i)t 

{Og[t) Og[t + T)) e ^^^^^^ _ i) + 4(2iV + 1){2N + 3) 

I N{N-l) _ 
4(2iV- l)(2iV+ 1) 4(2iV+l)2 / -l - ^ 

These expressions yield for r = 0, t — > oo the CUE results 

N iV^ 

(^)CUE = y ' '\^9 )CUE = 4(4jy2 _ (2.51) 
and for r = t = the COE expressions 

(^7)coE-y-2(2Ar+l) ' ('^^ ^^^^ " (2iV + l)2(2iV + 3) ^^'^^^ 



which are given in Ref. ||11|| . 

In the limit oo, one has to keep constant the ratio t/tc = 8Nt in order to get a non 
trivial crossover. Eqs. ( |2.49D and ( ^.5U| ) then become up to corrections of order 1/A^ 

1 

{9{t)) = J - ^e-'/'^ , (2.53) 
{Sg{t) 5g{t + r)) = ^ fi-^/*^ (l + e^^*/*.) . (2.54) 

The average ( p.53| ) can also be obtained by the asymptotic density ( |2.45| ). 

Fig. 5 and Fig. 6 give illustrations of the expressions ( p.49|) and ( p.50|) for the average 
conductance and its variance at = 1, 5, 00. 



III. 5-BROWNIAN MOTION TIME AND MAGNETIC FLUX. 



Till now, we have considered the Brownian motion ensemble COE CUE in terms of 
the parameter t. We suppose as in Refs. ||10|JTT[] that the scattering matrix of a chaotic 
quantum dot is well described by a COE or a CUE matrix for zero or sufficiently large 
magnetic flux respectively. Then the obvious question concerning the relation between the 
Brownian motion time t and the magnetic flux $ arises. This section is devoted to this 
question. The consideration of the magnetic flux requires some microscopic model. Usually, 
such models are treated by numerical or semiclassical means P J2^J23[] . Here we use another 
approach that connects S with a Pandey-Mehta Hamiltonian describing the crossover 
between gaussian ensembles for the energy levels (COE CUE). The crossover parameter 



of the latter has been related to the magnetic flux for disordered systems [21] and for ballistic 



cavities 23,22 
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The argumentation of Ref. p2| rehes essentially on numerical calculations and that of 



Ref. 



23[ on a semiclassical approximation. In principle, the relation between the Pandey- 

However, 



Mehta parameter and the magnetic flux is understood due to these references, 
we reconsider this relation in subsection |1II A| , using arguments that account for both the 
disordered and the ballistic case. The known result for the disordered case of Ref. ||21[ 
is recovered. In the ballistic chaotic case, we find that the relation between the Pandey- 
Mehta parameter and the magnetic flux is also governed by the quotient of two characteristic 
energies. In subsection pi B| we discuss briefly the typical field scales of a ballistic billiard. 
In particular, for a very strong magnetic field the system becomes integrable due to the 
appearance of Landau-levels. Finally, in subsection [111 C| , we relate the Pandey-Mehta model 
for the Hamiltonian with the S-Brownian motion ensemble for the scattering matrix. We 
find that, in order to map the S'-Brownian motion model onto this more microscopic model, 
we need to assume a certain relation between the Brownian motion time and the Pandey- 
Mehta parameter. This identification gives identical expression for the average quantities 
in the limit of small fluxes and large number of channels, but does not solve remaining 
discrepancies for the correlation functions. 



A. Magnetic flux and Pandey-Mehta Hamiltonian 

The generic random matrix ensemble that describes the crossover from GOE —>■ GUE is 
the Pandey-Mehta Hamiltonian 

Ho, = Hgoe{v^) + iaA{v^) (3.1) 

where Hqqe is a GOE matrix {v^ being the variance of the non diagonal elements) and A 
is an antisymmetric matrix whose independent matrix elements Aij are gaussian distributed 
with variance f ^. a G [0, 1] is the crossover parameter, a = corresponds to the GOE-case 
whereas a = 1 is the GUE-case. 

A typical one-particle Hamiltonian is of the form 

H = — (p-elY + V{f) (3.2) 
2m ^ ^ 

in d = 2 or (i = 3 space dimensions where A = ^{B x r) is the vector potential in the 
symmetric gauge. The potential V{f) can model the boundary of a chaotic ballistic cavity 
or the impurities of a disordered conductor with bulk diffusion, such that the energy levels 
have Wigner-Dyson GOE-statistics without magnetic field. We decompose if in a real and 
imaginary part, H = Hi+iH2 with = §^ + V{f) + j^A\ H2 = -^p- A = -^B ■ L und 
L = f X p is the angular momentum. The contribution ~ A"^ in Hi modifies the potential. 
For the purpose of this subsection we restrict ourselves to the case where the magnetic field 
is sufficiently small and does not change the GOE statistics for Hi ^ H{B = 0). The range 
of magnetic field at which this behavior is not true corresponds to a strong influence of 
the Landau Levels and the system can not be considered as "chaotic". A more detailed 
discussion of this topic and the relevant scales of the magnetic field will be given in the next 
subsection. 
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We consider a certain energy Ep and want to describe H by the Pandey-Mehta Hamilto- 
nian (|3.1| ) in an energy range E G [Ep—AE, Ep+AE]. When B = Owe assume that H = Hi 
with HQOEiv'^)- The variance is determined by the requirement v = AqVCV" /it where 



Aq (-C AE) is the level spacing of Hi and A/" is the chosen dimension of the GOE ensemble, 
assumed to be large. We do not discuss here whether the Pandey-Mehta Hamiltonian gives 
an appropriate description of the microscopic Hamiltonian H a.t B ^ 0. The numerical and 
perturbative calculations of Ref. |2T[] confirm that the level correlations of a disordered ring 



in the diffusive region agree with the correlation functions of Ha found in . The parameter 
a is then given by the identification 



-=J^r^(2^lr] (3.3) 



AT Ao V $0 



where Et = TiD / {2tiRY {D = diffusion constant, R = ring radius) is the Thouless energy 
and = h / e is the flux quantum. The dependence on A/" indicates that the GOE to GUE 
crossover is properly described by the parameter A = VJ7 a/vr in the limit A'' — oo . The 



ballistic quantum billiard has been studied in Ref. ||22| by numerical diagonalization. Again 
the correlation functions of are valid and the parameter a is found to be proportional to 
a flux $ applied through a typical of the system. 

We now give an alternative derivation of (|3.3| ) based on a crude estimation of the variance 
of a typical matrix element {H2)ij in the eigenvector basis of the real part Hi of H . Let \i > 
be an eigenvector of Hi with eigenvalue For Ei ~ Ep, we consider the expectation value 

< i\L,{0) L,{t)\t >=< i\L,e'"^'/^L,e-'^'*'^\i >= ^ |< t\L,\j >f e*^^^-"')*/^ (3.4) 

j 

where is the 2;-component of the angular momentum (assuming that the direction of B 
is chosen as the 2;-axis). Lz{t) is the time dependent angular momentum with a dynamic 
determined by Hi. We now average over the disorder or over the random shape of the cavity 
and assume that the average of the matrix element of the angular momentum is a smooth 
function of the energy difference Ej — Ei, i.e. <^|< i\Lz\j >|^)^^ — K^j ~ ^i)- Eq. ( p.4|) then 
becomes 

{LMLMq,v = ^Jd^ ^(^^) ^2 {^) e^-* . (3.5) 

The subscript on the left side denotes a combined expectation value from quantum mechanics 
and from the random parameters of V. R2 is just the two point correlation function in the 
orthogonal case with the limit R2{x) ~ 1 if ^ 1. The inversion of the Fourier transform 
gives 

lihu) R, (^) = ^J dt {LMLMgy e'^^* • (3.6) 

The expectation value in the integral vanishes at a time scale \t\ ^ where is the corre- 
lation time for the angular momentum. Hence, the integral itself vanishes as a function of lj 
at the scale \uj\ ^ 1/tl. From this behavior we can already see the limitations of the GOE 
description for Hi. The absolute squared matrix elements of an arbitrary fixed operator with 
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H 



respect to the eigenvectors of a GOE-Hamiltonian yields always a constant value indepen- 
dent of the particular indices i ^ j after the average, since there is no preferential basis. Eq. 
( p.6| ) implies that for Hi this is only true for an energy difference \ej — ei\ -C H/tl =■ Ecow 
This limits the size of the energy range that we may consider. Since we are interested by 
the off diagonal matrix elements for Aq ^ \ej — Si] <^ -Ecorr, we take i?2(^) — 1- From Eq. 
]6D we get for the variance of the imaginary matrix elements (-^2)4^ = < ^I-^zIj > of 

where the correlation time is defined as {Ll)Qy = ^ J dt {Lz{0)L;,{t))Qy. Let 2R be 
the typical diameter of the system and fi;' be the Fermi velocity. Furthermore, we introduce 
a characteristic ballistic energy scale -Ebaii = fivp/i^R) and the flux $ = ttR^B applied 
through a circle of radius R. The expectation value for the angular momentum is just 
{LDqy = Km^VpR'^ where k is a numerical factor characteristic of the shape of the dot 
which can be evaluated in a semiclassical way, as discussed in Appendix y. 
The variance (|3.7|) becomes now 



^corr 



and has to be compared with the variance of the imaginary matrix elements in ( |3.1|) that is 
just a^AgA/'/vr^. This gives for the parameter a 



. (3.9) 



AT VAo ^corr V$o 

This expression contains the dimensionless flux and three typical energies -Ebaii = hvp/i'iR), 
-f^corr = ^/tl and the level spacing Aq. 

Let us first consider a disordered conductor where the electrons are scatterered by many 
impurities. The angular momentum correlation time is then the elastic scattering time 
tl = Tel and the energy quotient becomes i?baii/(-^corr Aq) = {2'kY Et d / /S.o where Et is the 
Thouless energy {Et = TiD /{2ttRY for a diffusion constant D = VpTe\/d). Eq. ( p.9| ) confirms 
the result (|3.3| ) of Ref. |^ for a ring with k = 1/d (Appendix ICj). 



For a ballistic quantum billiard, equals the time of flight k2R/vf {k is a further 
numerical factor) through the cavity and -Ebaii = ^Ecow This gives 



^ 4^^/* . (3.10) 



V N Ao V$o 

Both, this result and Eq. ( p. 3D for the disordered case depend on the energy ratio Eq/Aq 
where tc = h/Ec is the typical time for the electron to cross the system, either by a direct 
ballistic or by a diffusive motion. 

It is well known [^] that the level statistics of a disordered conductor are well described 
by the GOE-statistics up to an energy scale Ep- Our derivation yields a smaller energy scale 
-E'corr = ^/tci- This is due to the fact that we use a criterion based on the matrix elements of 



15 



the angular momentum. If one estimates in a similar way the matrix elements of the position 
operator, one recovers the characteristic scale Ex-, which matters if we consider only for the 
level correlations.^ 

Such a difference does not matter for a chaotic ballistic dot where those two energies 
merge to a single characteristic scale E'baii = hvp/i^R). The energy ratio -Ebaii/Ao is of the 
order of [kpRY^^ ^ 1 {kp is the Fermi momentum, for comparison: Ep/Ao ~ {kpRY)- 
The critical flux for the GOE to GUE crossover (for an energy interval of order Aq ) is 

*c ~ ^o/ikpR)^'^-^^'^ < $0- In ci = 2 dimensions one gets Je^~^o ~ {kpRf^ ~ n^'^ 



where n is the number of energy levels below the Fermi energy in the quantum billiard. 
This scaling behavior was already given in Ref . p2| and applied in Ref . [§] . In the numerical 



simulations of Ref. the factor y -Ebaii/ is of the order of one because it is rather difficult 
to increase n to values where this scaling behavior becomes important. The situation changes 



a bit in three dimensions, where yi?baii/Ao ~ (kpR) ~ n^^^. An estimation of the numerical 
factor that appears in Eq. (|3.10|) is given in Appendix 0. 



B. Strong and weak magnetic field scales in a quantum billiard 

We now briefly discuss the effects of stronger magnetic flelds in a quantum billiard. 
For simplicity we consider the two dimensional case, where a free electron in a magnetic 
field B has Landau levels with energies En = TioJciji + \). n = 0, 1,2, . . . is the number 
of the Landau band and uoc = {eB)/m is the cyclotron frequency. The eigenstates are a 
combination of harmonic oscillator states (with frequency Uc) in one direction and plane 
waves in the orthogonal direction. The magnetic length Ic = \Jfi/ {mujc) is the characteristic 
length of the harmonic oscillator. Ic is the typical width of the states (in the oscillator 
direction) of the first Landau band with = 0. At an arbitrary Energy Ep the states have 
a width ~ ^ Ep / (huc) Ic ^ vp/uc = R\am of the order of the Lamor radius. 

If in a quantum billiard of typical size R the magnetic field is very high so that -Riam R, 
the effect of the Landau Levels becomes very strong. There is a finite number of states of 
each Landau band (with energies < Ep) that do not feel the boundary because their typical 
width -Riam is sufficiently small. The states near the boundary ( "edge states" ) are raised to 
higher energies. The energy spectrum near Ep does not differ very much from the spectrum 
of the unbounded plane and cannot be described by a GUE spectrum. The only effects are a 
(very) small breaking of the degeneracies and the appearance of the edge states. The classical 
trajectories are circles with radius -Riam that do not see the boundary. If the magnetic field 
is lowered (such that Ic ^ R ^ -Riam) the fast electrons at the Fermi energy Ep feel very 
well the boundary and the energy spectrum should be closer to a GUE statistics. But at 
low energies E < Ep the Landau levels still exist and the corresponding edge states can 



^ The considerations made here for the disordered case are only valid in the diffusive or metallic 
regime where the states are not localized. In the localized regime the GOE statistics are no longer 
valid and, in addition, the expectation value of would depend on the localization length instead 
of the typical system size. 
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have energies near Ep. The latter may well complicate a simple GUE picture because the 
corresponding eigenfunctions have a very small overlap with most of the states at Ep. This 
effect disappears if the magnetic field becomes smaller such that R < Ic- Then the real GUE 
statistics apply. 

Taking into account the effects described above and the results of the last subsection, we 
distinguish four regions for the magnetic field. 



(i) $/$o < yAoZ-Ebaii, the GOE to GUE crossover regime. The level statistic of this 

regime is described by the correlation functions of Ref. 0. If $/$o ^ ^Ao/-Ebaii, the 
behavior is near to a GOE but one must take into account that the level repulsion 
becomes immediately quadratic on a very short energy scale. 

(ii) y/AoZ-Ebaii ^/^o ^ 1) the GUE regime. The level statistics are described by a GUE. 

(iii) 1 <^ ^/^o <^ Ep/E\,si\ or Ic ^ R ^ -Riam, the regime where first Landau level effects 
appear. The level statistic is roughly GUE with complications due to edge states. 

(iv) Ep/E]~,ai\ <^ $/$0) or -Riam ^ -R, the Landau level regime. The behavior of the 
eigenvalues near Ep is very close to the spectrum of the Landau levels and thus regular. 

This subdivision of the different magnetic field scales remains essentially valid in three 
instead of two space dimensions. The behavior of the regions (i) and (ii) is the same. In the 
regions (iii) and (iv) the influence of the ^-component of the momentum becomes important. 
In the region (iv) every energy eigenvalue is a sum of one Landau level of the xy-plane and a 
contribution due to the kinetic energy in the z-direction. The Hamiltonian is (roughly) a sum 
of two regular independent parts. The situation in region (iii) is again rather complicated 
due to a mixture of "Poisson"- and "GUE" -levels. 

Our main concern is the crossover between the regions (i) and (ii) where we suppose that 
the Hamiltonian is well described by the Pandey-Mehta Hamiltonian (|3.1|) with a parameter 
a given by Eq. ( |3.10| ). In the next subsection we discuss how a can be related with the 
Brownian motion time t of the GOE GUE crossover. 



C. S'-Brownian motion time and Pandey-Mehta Hamiltonian 

The 2N X 2A^-scattering matrix S of a physical system that is described by an A/" x M- 
random Hamiltonian H can be expressed by 



S{H) = l-2m G+ W 
where we have used the abbreviation 

G± = {E - H ±t7TWW^)-^ 



(3.11) 



(3.12) 



W is a N" X 2N matrix that describes the coupling between the states of the system with 
the scattering channels. The expression ( 3.11 ) is well justified (and references therein) 

and is the usual starting point for supersymmetric calculations |17|,|8| J28|J29[| . It is not very 
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difficult to see that S{H) given by is indeed a unitary matrix if H is hermitian. We 



take the hmit A/" — oo but the dimension 2N of the scattering matrix may remain finite (as 
considered in Ref . [§] ) . When we consider the hmit ^ 1 , we nevertheless assume ^ Af. 
We note M = 2N the dimension of 5*. 

The expression ( |3.11| ), relating an appropriate Hamiltonian to S, is what we need for 



giving a microscopic motivation of the S'-Brownian motion model. However, we assume for 
simplicity a crossover random matrix model for the Hamiltonian, and not a real microscopic 
model. 

Let us assume that the Hamiltonian H contains a typical parameter depending on the 
magnetic field. A small change of the parameter leads to the replacement H H + AH 
which gives for S 

S{H + AH) = S{H) exp(-z5X) (3.13) 



where 6X is a hermitian matrix. Expanding Eq. ( |3.11| ) up to second order we obtain after 
some matrix-algebra the expression 

6X = 2tt[w^ G_ ah G+ W 

+W^G-AHG+{E-H)G-AHG+w} + 0[{AHf) . (3.14) 

We now assume that the Hamiltonian can be written in the form 

H{^) = Ho + ^Hgve{v^) (3.15) 

where Hq is some (fixed or random) part of the Hamiltonian that may be rather general 
at the moment. -ffGUE('y^) is a GUE-random matrix {statistically independent of i^o) whose 
elements are independent and normally distributed. The average and variance of Hque{v'^) 
are characterized by 

{Hgve{v^)) = , {HcvEivyAHcMv^)) = tT{A) 1^ (3.16) 

where A is an arbitrary JV x J\f test-matrix. We consider now an infinitesimal change of the 
parameter 7 — > 7 + ^7. Since the sum of two normally distributed variables is also normally 
distributed, we can apply the decomposition 

{-f + Hgve{v^) J HgveAv^) + A-f HgveAv^) (3.17) 

where i^cuE,! and i^cuE,! are two different and completely independent GUE-matrices both 
characterized by ( |3.16| ). The variances are related through 



(7 + 57)2 = y + (A7)2 ^ A7 ~ ^27 57 . (3.18) 

The replacement ( |3.17|) is exact and does not incorporate any approximation if the trans- 
lation ( |3.18| ) is properly used. We emphasize that this substitution of one set of Gaussian 
variables with a sum of two independent sets of new Gaussian variables works as fas as 
we are only interested in simple but arbitrary averages of the type {f[S{H{'y))]) (where 
/ is an arbitrary function of the scattering matrix) for a given value of 7. In the case of 
correlations between different 7 values, i. e. if we want to calculate averages of the type 
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(/i[S'(if(7i))] /2[5'(-ff(7i + 72))]), the decomposition ( p.l7|) can not be used because it is not 



possible to vary 72 in this way without changing 71. However, for the case 71 = the de- 
pendence on 72 can be considered in the same way. Therefore, the following argumentation 
concerning the relation between the Brownian motion time t and the parameter 7 holds for 
averages {f[S{H{'y))]) and the particular correlations {fi[S{H{0))] f2[S{H{'y))]) . 

We identify the small matrix AH with A7iJGUE,2- The hermitian matrix 6X given 
by Eq. ( |3.14D is now a random matrix matrix because of i^GUE,2- We average over this 
GUE-matrix (denoted by (■ ■ ■)2) and obtain with the help of ( p.l6| ) 



(5X)2 = ^i;2(A7)^(tr(G+)+tr(G_))Q(iJ) , (3.19) 
{6X^ A6X)2 = ( A7)2 tr(i Q{H)) Q{H) (3.20) 

where now A is an arbitrary M x M-test matrix and Q{H) is the Wigner time delay matrix 
given by 

(99 

Q{H) = 27iW^G.G+W = -iS^ — . (3.21) 

oE 

We underline that the matrix elements Qij give the time delay characterizing the scattering 
modes of the chaotic cavity. This physical meaning is very useful for a more concrete and 
intuitive understanding of this somewhat abstract S-Brownian motion ensemble. We have 
to compare the averages ( p.l9| ) and ( p. 201 ) with the corresponding averages of the Brownian 



motion "increment" SX^^^ |15| , p!9| given by 

(5X(^))2 = , (3.22) 
((5X(^))t A 5X(^))2 = D5t tr(i) Ui (3.23) 

where 6t is the change of the Brownian motion time t and D is the diffusion constant (which 
is chosen as D = 4 in the applications of this paper and in Ref. [Il9|). From this we can 
directly state: 

The Brownian motion model is the idealization in which the following two conditions are 
fulfilled 

tr(G+) +tr(G_) = , (3.24) 
Q{H) = go(7') 1m (3.25) 

where qoi'j^) is a number that may depend on 7^ but not on the different realizations of 
Hque (|3.15| ). If these conditions are valid the relation between the Brownian motion time 
and the parameter'-) is determined by D 5t = f ^(A7)^ 50(7^) ~ 27^7 55(7^) which results 
in the translation 

t = ^ r dx qlix) . (3.26) 
D Jo 

Assuming ( p.25| ) we can evaluate go (7^) as the trace of the matrix Q{H) which is just 
Wigner 's time delay 
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M ■ go(7^) = tT{Q{H)) = (-i)tr(G+(27ri WW^)G-) 

= (-0(tr(G_)-tr(G+)) (3.27) 



where the last equahty holds due to ( p. 12 ). Let us assume (as in Ref. [^]) that the operator 



H + mWW"^ has complex eigenvalues + zS^, > such that 

The second equality is an approximation valid if the typical values of the imaginary parts 
satisfy 3> A where A is the level spacing of H . Then, one can replace the /x-sum by an 
integral over E'^. This condition is related to the limit M ^ 1 where the conditions (p.24|) 
and (|3.25| ) make sense. Actually, Eq. ( p. 2^ ) yields already the first condition (|3.24|) . 



The second condition ( p.25| ) means that the time delays of the different scattering modes 
are roughly identical (which looks natural for a chaotic cavity) and do not fluctuate (in 
comparison with their mean value). This second assumption is less natural, but can be 
justified if the number of channels is large.0 

In the following, we only consider a range for the parameter 7 that is small enough for 



having A determined by the part Hq of the Hamiltonian ( p.l5|) , i.e. A Aq where Aq 
is the level spacing of Hq. This assumption is very reasonable since in the large M limit 
the range 7 ~ l/vCV" yields already the wanted crossover (if Hq is a GOE matrix). The 



number go is now approximately independent of 7. Due to (|3.28|) and (|3.27|) it is given by 



go — 27r/(MAo). Furthermore, since the variance is related to the level spacing Ague of 
-f^GUE via ||2^ = A/" Aque/tt^, one gets for the Brownian motion time (|3.26| ) 



■J^ ^GUE 2 /o2q^ 

M2 Ag ^ ^'^■^^> 



where the diffusion constant D = A. Assuming the validity of the conditions ( fj.24D and 



( p.25| ) with the approximation (|3.28|) , the S-matrix of H = Hq + •y HQUEiv"^) is described 
by the Brownian motion ensemble with an initial condition S{t = 0) = S{Ho) and a time 
parameter ( |3.29| ). 



Now, we assume that the initial Hamiltonian is a Gaussian orthogonal random matrix 
Hq = Hgoe{v'^) and that H{'y) = Hqoe{v'^)+'J Hgue{v'^). The S'-Brownian motion ensemble 
corresponds approximately to this choice with the S'-matrix expression (|3.11|) and the time- 
translation if: = (A/'7^)/M^. The Brownian motion ensemble can describe properly averages 



^ If Hq is a GOE random matrix, we can consider the deviation 5Q = Q — 1/Mtr((5) 1m from Q 
with the matrix qqIm- A measure for the fluctuations 6Q is just tT{6Q^) = tr{Q^) - l/Mti{Qf. 
We have calculated the average of this quantity using supersymmetric technics in leading order 
of a large channel number with the result (tr((5Q^)) = (27r/Ao)^/M = 0{M~^) where Aq is the 
level spacing of i/cuE- We do not present this calculation based on a perturbative treatment of 
the cT-model [^8|. The relative variance {\SQij\'^) / {Qu)'^ of an arbitrary matrix element of 6Qij is 
of the order 0{M~^) thus verifying the second condition ( |3.25 ) in the limit M ^ 1. 
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and correlations between 7 = and 7 > 0.0 

We can note that the random Hamiltonian H{'y) is not identical to the Pandey-Mehta 
Hamiltonian However, we can again use a decomposition of Gaussian variables which 

leads to the identification — > Vl — «^ Hqo^{v'^) + \/2a Hq\]-e{v'^). The factor a/2 is due 
to the fact that the sum of the real and the imaginary parts of the (off diagonal) matrix 
elements has a doubled variance, i.e. ((ifGOE.jj)^ + {.^ijY) = = 2(|ifGUE,jjP)- Since 
the typical scale of the parameter a is a ~ 1/VJ7, we can approximate Ha — HQOEiv'^) + 
\/2 a HGUEiv"^) = H{y/2a) and we obtain thus the relation 7 = \/2a leading to the final 
relation between t and a 

^ = ^(«) = ^(A^"') = ^-'A^ (3-30) 

where A is the crossover parameter used in Ref. ||^. 

We emphasize that the identification of the Pandey-Mehta Hamiltonian with H{\/2a) 
works only as fas as averages are concerned. The correlations between a = and a > 
are then not properly described by the Brownian motion model and indeed the critical field 
scale for these correlations obtained from the Brownian motion ensemble is a factor 1/2 
too small if compared with semiclassical results and a supersymmetric treatment of the 
Pandey-Mehta Hamiltonian [PD|. 



IV. VALIDITY AND LIMIT OF THE 5-BROWNIAN MOTION ENSEMBLE 



The S'-Brownian motion time t is related to the crossover parameter A = ^/Maj-n of the 
Pandey-Mehta Hamiltonian ( |3.1| ) 



via 



t(A) 



27r^A^ 
M2 ■ 



(4.1) 



This identification assumes the limit M ^ 1. One can express A in terms of the magnetic 



flux: for a ballistic chaotic dot, we found in subsection |III A| the expression 



/ 4k k hvp ( $ 



(4.2) 



^ We have still to ensure that the initial condition 7 = corresponds to the circular orthogonal 
ensemble. In general, the S-matrix described by Eq. ( 3.11| ) represents a wide class of random 
matrix ensembles depending on the independent matrix elements of the symmetric M x M-matrix 
W'tVF §. In the case H = Hqqe^ one expects that there is one particular choice of W^W in 
which S{H) is exactly distributed as a unitary COE random matrix. This choice is characterized 
by the requirement that the average {S{H)) = vanishes corresponding to a maximal coupling 
between the scattering channels and the system ||7). At the energy range l-Bl <C AAAqoe one 
requires ||^J2^ W^W = MAqqe/'^'^ • 1m- The typical value of the imaginary parts used in Eq. 
(|3.2§| ) can then be estimated as ~ 1/M tr(7rVFVK^) = MAqqe/t^ and the condition for the 
validity of (|]2|) reads as M > 1. 
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where k and R are numerical constants characteristic for the shape of a dot of typical size 
2R and of level spacing Aq without magnetic flux, f is the Fermi velocity. 

We can now compare the result ( p.49| ) for the average magneto conductance with that 



of Ref. P). From Eq. ( |4.1|) we obtain for the S'-Brownian motion ensemble 
whereas the supersymmetric result of Ref. P] reads 

Eq. (^4.4|) is not the exact analytical result, but an accurate fit obtained after the numerical 
evaluation of three integrals In leading order of 1/M and in the region -C M/(87r^), 
both results become identical whereas for larger values of A the difference between the 
Lorentzian and the Gaussian is more important. This can be understood from the argumen- 
tation developed in subsection ^11 C| , which shows that the S-Brownian motion model cannot 
be exactly mapped on the ^-matrix expression ( |3.11D used in Ref. . The expressions ( |3.19|) 



and ( p. 20 ) for the averages of the increment 5X show that in reality we have to deal with 



a more complicated stochastic process than the simplified one implied by Eqs. ( p.22| ) and 
( p.23| ). Both stochastic processes are only approximately in agreement for sufficiently small 
Brownian motion times or values of the symmetry breaking parameter, i.e. for A^ <^ M/%, 
where the results ( [4.3[ ) and ( [4.4[ ) are consistent. 

The correlations between two different values A and A + AA are only properly described 
if A vanishes and AA is sufficiently small. In addition, the comparison with semiclassical 
and supersymmetric [jl9[ results gives a factor \/2 for the relevant AA scale. This factor 



can be understood by the difference between a symmetry breaking due to an imaginary 
antisymmetric matrix or to a fully hermitian matrix (compare subsection [111 C| and Ref. 

In Ref. [|T^ it was shown that the Fokker-Planck equation for the transmission eigenvalues 
Tj can also be obtained for a more general matrix 6X than described by Eqs. (|3.22| ) and 
( p.23| ). If we assume that this is true for the S'-matrix ( |3.11| ) with the Pandey-Mehta 



Hamiltonian, we have only to determine correctly the relation between the S'-Brownian 
motion time and the parameter A. Then, the expression ( p.30[ ) will only be an approximation 
valid in the range A^ <^ M/ (87r^) and one could use the supersymmetric result ([4.4|) to obtain 
a more appropriate expression 

1 , / 8(M-l)7r2An 

In 1 + ^ ^ • (4-5) 



AM \ M2 ^ 

This modified relation assumes that the problem resulting from the difference between the 
S'-Brownian motion model and the microscopic S-matrix ( |3.11| ) can be solved by an appro- 
priate rescaling of the Brownian motion time.0 This will work for the average conductance 



^ In fact, the relation ( [4.5| ) gives the correct A-dependence of the conductance fluctuations (in the 
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(suppression of the weak localization correction), but not for the conductance- conductance 
correlations. 

It is worth to underline that the crossover from orthogonal to unitary symmetry happens 



at the critical M- dependent value Ac = yM/S /it for the parameter A. The corresponding 
critical flux is given by 



M Ao2R , , 

<!>c = <l>o J -—^ ^— . 4.6 

This is due to the M-dependence of the typical time for an electron to stay in the dot. With 
increasing channel number N = M/2, this typical time is shorter. Therefore the critical flux 
needed to break the phase coherence between time reversed paths is larger. The semiclassical 
analysis of the Ref. |Q shows that the distribution of the effective area 6 surrounded by the 
electronic trajectories plays an important role. The comparison (of the quadratic term) of 
the Lorentzian of p with the above expressions ( |4.3| - |4.4D yields just 



<l>c = ^a,iA% (4.7) 

where A is the area of the two-dimensional cavity and ad is the characteristic inverse area 
appearing in the exponential behavior of the probability density p(0) ~ exp(— ctd \ for 
e. Eqs. (U) and (j^) yield directly ad ~ or (|e|) ~ 1/Vm. 

The studied S-Brownian motion ensemble has the decisive advantage that many impor- 
tant properties (Section JJ) can be calculated exactly, even for a finite dimension M of the 
S-matrix. It is possible to relate this model to the Pandey-Mehta Hamiltonian when there 
is a large number of channels up to a certain approximation and to express the Brownian 
motion time in terms of the parameter A ( directly related to the magnetic flux). This 
yields consistent results for the magneto conductance for sufficiently small values of A -C A^ 
whereas in the other limit the deviation between a Lorentzian and a Gaussian behavior 
becomes relevant. This indicates that the S"- matrix ( 3.11| ) with the Pandey-Mehta Hamil- 



tonian ( |3.1| ) is not equivalent to the S'-Brownian motion model. This latter model assumes 
a uniform and irreversible diffusion process in the S'-matrix space being valid for short time 
scales. For longer times, a rescaling of the time scale, as proposed in Eq. ( [4.5|) , improves 
the validity of the model, as far as averages are concerned e. g. average and variance of 



the conductance. In Ref. |3^, it is however shown that even this rescaling of the Brownian 
motion time is not sufficient for having the right autocorrelation function between different 
times. 

To improve this exactly solvable model, one might reconsider its two main assumptions: 
first that the infinitesimal evolution exp{—i6X) of S is statistically independent of S and 
second that this evolution is isotropic with a time delay matrix Q{H) fulfilling the conditions 
Eq. (p.25|) and ( p.24|) . In principle, these conditions are only valid in average. To remove the 
first hypothesis would imply a Brownian motion with memory effect or a diffusion constant 



limit M ^ 1) if we apply this relation to ( 2.54 ) (for r = 0) and compare it with a perturbative 
treatment of the cr-model . 
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becoming a random quantity on itself: i. e. a rather involved Brownian motion in a random 
medium. 

It might be interesting for future works to use the S-Brownian motion ensemble for 
a probably simpler problem: the parametric correlations with respect to the energy. In 
this case, the time delay matrix Q{H) appears in a very direct way as 5X/5E and just 
the fluctuations around its average (which we have neglected in the last subsection) are 
decisive to generate the Brownian motion dynamics. One could numerically investigate if 
the statistic of Q{H) in a chaotic cavity is at least approximately independent of the statistic 
of S and obeys the isotropy condition necessary for our model. Furthermore, one could try 
to formulate a non-isotropic Brownian motion if imposed by a more realistic time delay 
matrix. 
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APPENDIX A: CALCULATION OF THE SKEW ORTHOGONAL 

POLYNOMIALS 

In this Appendix, we sketch the calculation of the skew orthogonal polynomials q'^\x) 
for the antisymmetric scalar product (|2.11|) . First we consider the case t = 0. Due to 
g{x, y; 0) = 6{x — y) and ( |2.7| ) the scalar product then becomes 

<f,9>R :=< /, 9 >^°)= l C dy, [' dy,^h=^^= f{y{) g{y,) . (Al) 

In addition, we need the conventional symmetric scalar product given by 

<f.g>c= l\f{y)g{y) (A2) 

with the Legendre polynomials Pn{x) as orthogonal polynomials, i.e. < Pm >c= 
2/(1 + 2n) Snm- In the following we omit the superscript (■ ■ •)*'°^ ^-nd denote qn{x) the 



skew orthogonal polynomials (compare Eq. ( |2.12|) ). In usual random matrix applications 



the two type of scalar products contain the same weight function and one can exploit 
a simple symmetry due to parity x —x. The skew orthogonal polynomials are then di- 
rectly related with the conventional orthogonal polynomials and their derivatives [^. The 



situation here is a bit more complicated and our aim is to expand the qn{x) in terms of the 
Legendre poljTiomials. We use therefore the ansatz 

n 

qn{x) = X! ^nmi(l + 2"^) Pm i^X 

) (A3) 

m=0 



24 



where the bnm are coefficients to be determined. We note that both scalar products are 
related via the following Lemma 

< f,g(x)-gi-l)^^^ >^ = ^ < f,gix)-2il-x)g'ix) >n (A4) 



where x is the argument of the function on the right sides. The Lemma ( |A4| ) is easily verified 
by two times partial integration of < f,g (and using that e{u) = sign(M)). Using the 
expansion 

1 oo 



v/2(l-x) t 

and standard identities which relate the Legendre polynomials and their derivates we find 
for an arbitrary function f{x) and for g{x) = Pm{x) from (|A4D 



<f,p^-(-irY.Pk>c= 

k=0 

= ^ < /, (1 + 2m)Pm - 2{2m - l)Pm-i + 2P'^^^ + 2P'^^^ >r . (A6) 

The coefficients bnm can be evaluated by bnm =< Qn, Pm >c which follows directly from ([A3D. 
This identity holds also for < m if we put in this case bnm = 0. The skew orthogonality 
of the Qn implies < Qny^m >R= if is an arbitrary polynominal of degree m < n — 2 (or 
m < n if n is an even number). We put f = qn if n > m and obtain from (|A6|) , (|A3|) 

n ^ 

bnm - i-ir T.^nk = -< Qn, (1 + 2m)Pm " 2(2m - l)P„_i >R (A7) 
ifc=0 ^ 

where we have used that -P^-i Pm~2 polynomials with a degree <m — 2<n — 2. 
In general this equation holds for m < n but in the case of even n it is valid for m = n + 1, 
too. 

We consider ffist the case of even n = 21. If m < 21 the r.h.s. of ( [AT] ) vanishes and we 
get therefore 

21 

b2i,m = {-ir E b2i,k = {-irb2i,2i (A8) 

k=0 

where the second equality holds because the sum does not depend on m. For m = 2/ + 1 
Eq. (^ gives 

21 2 

b2l,2l = E ^2Z,fc = IT < q2h (3 + 4/)P2«+l >R= (&2«+l,2/+l) • (A9) 
A:=0 ^ 

The last equality is valid because the Legendre polynomials can be expanded in terms of 
the Qn by inversion of ( [A3|) . The skew orthogonality implies that for ( [X9| ) only the ffist 
coefficient of this expansion is needed. 



25 



We exploit now (^3) ^^^^ ^"^"^ n = 2/ + 1 which is a bit more involved since 

three subcases m < 21, m = 21, m = 21 + 1 have to be considered yielding the following 
equations 

= (-l)"s2;+i , m<2l , (AlO) 

b2l+l,2l = S21+1 - 7 (All) 

021,21 

U I ^21+1,21 ,2 /A ION 

02/+l,2«+l = -S2«+l + 7 7 ^ 7 (A12) 

02Z+1,2«+1 02/, 2« 02/, 2/ 

where we have used the abbreviation 

2/+1 

S2/+1 = hi+i^k ■ (A13) 

fc=0 

For ( |A12|) the first two coefficients in the expansion of P21+1 in terms of the were needed. 
Eqs. (|A9|) , (|A11| ) and ( |A12| ) yield two independent equations for the four quantities 



&2/,2i, b2i+i,2h b2i +i,2i+ i and S2/+i- The other 62;,™, &2Z+i,m (with m < 21) are determined by 
Eqs. (^^) and ( |A1(]| ). We have therefore two free parameters that are related to the simple 
transformation q2i ^ A- q2i and q2i+i —>■ ■ q2i+i + B ■ q2i (with arbitrary constants A 7^ 0, 
B) that does not change the skew orthogonality of the g„. Hence, we can choose 62Z,2« = 1 
and fe2Z+i,2i = which resuh in 62«+i,2«+i = 1, S21+1 = 1 and 62i,m = &2«+i,m = (-l)™ for 
m < 21. With these values of the coefficients b^m the ansatz ( |X^ ) yields the skew orthogonal 
polynomials for < >(*) at the time t = 0. 

It remains the task to find the skew orthogonal polynomials for an arbitrary time t > 0. 
We define now the polynomials 

n 

= E ^nm \{l + 2m) P^{x) e^"* (A14) 

m=0 



where Sn is given by = (1 + 2nY (compare Eq. (|2.4| )). The orthogonality of the Legendre 
polynomials and ( |2.4| ) imply directly that 

' dy g{x,y-t) q^\y) = q^^\x) = g„(x) (A15) 

-1 

with qn{x) given explicitly by ( |A^ ) and the above stated values of bnm- This relation shows 
immediately by virtue of Eq. ( |2.7] ) that the q^\x) are indeed a set of skew orthogonal 
polynomials with respect to ( p.ll|) . The result ( |A14| ) is just Eq. ( p.29|) stated in Section 

m 

APPENDIX B: MODIFICATIONS IN THE QUATERNION FORMULATION IN 

THE CASE OF ODD DIMENSION 



The quaternion formulation as described in Section [11 B| is only valid for the case of an 



even Number of channels. Here, we briefiy sketch the modifications for odd A^. It turns 
out that this case can be treated in a similar way as the even A^ case, but nevertheless 
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some technical complications appear. The following treating is rather parallel to the one of 
Ref. [|l7l, but adapted to the concrete situation given here. Since for an odd dimension the 
determinant of an antisymmetric matrix vanishes, Eq. ( |2.6| ) makes of course no sense if N 
is odd. The method of integration over alternate variables yields now for the integral 
the expression 



p{x, t) =oc \xi — Xj \ e 

i>j 



det 



-F{xj;t) 



1/2 



(Bl) 



where 



Fix;t) 



1 



dy 

1 ./2(l-y) 



g{x,y;t) = ^ P„(x) e 



-en t 



(B2) 



n=0 



g{x,y;t) is given by Eq. (|2.4|) and the second equality holds due to the expansion (|A5|) of 
the square root in terms of Legendre polynomials. H{xi, Xj] t) stands for an x subblock 
of an {N + 1) X (A^ + l)-matrix. F{xi] t) corresponds to the A^ + 1-th column and —F(^Xj] t) 
corresponds to the A^ + 1-th row. From Eq. ( |2.29| ) we get directly the identity 



dx q^^\x) F{x;t) = 1 . 

Therefore, it is useful to modify the skew orthogonal polynomials by 

q^Hx) = q^\x) - q^j^l.ix) , n<N-l 



(B3) 



(B4) 
(B5) 



In addition, we denote with Ql^\x) the dual functions that correspond to the q^\x), i.e. 
they are given by the integral ( p. 131) where the qn\x) are replaced by q^\x). Eqs. ( |B3[ - [B5D 
imply the "orthogonality" relation 



dx q^^\x) F{x-t) = 5n,N-i , n<N-l . 
Using this notations, we redefine the quaternion function a{x,y;t) by 

a{x,y;t) -- 



(B6) 



au{x,y;t) ai2{x,y;t) 
a2i{x,y;t) a22{x,y;t) 



(B7) 



with 



au{x,y;t) = -K-Q{x,y;t) + q^^^_^{x) F{y;t) , 

ai2ix,y;t) = Kgg{x,y;t) , 

a2i{x,y;t) = H{x,y;t) - KQQ{x,y;t) + 

+Q^S-MF{y-t)~F{x-t)Q^S-i{y) 
cr22{x,y;t) = KQ-{x,y;t) + F{x;t) q'ff_^{y) . 
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The function K is defined as in Eq. ( |2.15| ) but with the pecuharity that the sum runs over 
the range n, m = 0,1,. ..,A^ — 2. Eqs. ( 2.16 - 2.19 ) remain vahd, but in (|2.17|) one has to 
replace with — 1. Again the three main properties ( |2.21[ - pl23| ) are vahd and ensure the 
expression ( p.28|) for the m-point correlation function. The verification of (|2.21[ - ^l22|) is a 
rather lengthy but straightforward calculation using ( p.l6| - |2rT9| ) and ([B5|). The non-trivial 
property is again ( |2.23|) on which we will now concentrate. Instead of Eq. ( p.26|) , we find 
the identity 



Xi, Xj, t)) = limMi(5) ■ M2{S) 



<5->0 

where Mi^2{S) stand for the two 2A^ x 2A^-matrices 

'^'^ 6-Q^S_,{x,) + 6-'-F{xf,t) 



M,{5) 



q^XifZ 



( 



-5 ■ g(*) 



-Q{xj) 



N-l 



V H{xi,Xj;t) - 6 



[xj) + S 

' -2 



-'■F{xft) 
F{xi-t)F{xft) 



q[Xj) 




(B8) 



(B9) 



(BIO) 



with a quasi block notation similar as in Eq. (|2.26|) . The left column in Mi contains two 
A^ X (A^ — l)-matrices and the first row in M2 contains two (A^ — 1) x A^-matrices. The 
column (row) in the middle of Mi (M2) is of dimension 2A^ x 1 (or 1 x 2A^). The blocks in 
the right column of Mi (or in the lower row of M2) are just two square matrices of dimension 
N X N. The determinant of (IB 81) is now evaluated as 



det {A{a{xi,Xj;t))) = hm 

<5— »0 



Using the abbreviations 



(5^ ■ det (gi*^(xi))^ ■ det [-H{xi, xj; t) + ■ F{xi; t) F{xj; t) 



det 



a 




H{xi, Xj] t), Vi = 6 ^ F{xi] t) and the identities 
^ 



+ det 




V 



a 

T 



-V 




(Bll) 



(B12) 



(B13) 



we find that the expression in Eq. (pil|) become equal (apart from the normalization 
constant) to the square of the joint probability distribution ([Bl[). The second equality in 
( pi3|) holds because the determinant of an antisymmetric matrix in an odd-dimensional 
space vanishes. We have so far proved the last property ( |2.23| ), too. Hence, the m-point 
functions are again given by ( |2.28| ), but now with the quaternion function ([B7|). It turns 
now out that the results (|2]Mj2]4|) are valid for the odd A^-case, too. 



APPENDIX C: SEMICLASSICAL EVALUATION OF THE CONSTANTS k AND k 



In subsection [III A| , we introduced the numerical constants k and n that depend on the 



dot geometry. In this appendix, we give a semiclassical estimation of them for some simple 



28 



geometries. The constant k is defined through the relation {L^qy = nm^v'^R^ where the 
brackets denote a combined quantum mechanical and statistical average. In the disordered 
case the position r and the momentum p of the electron fluctuate on two different time 
scales tt ~ R^/D {D is the diffusion constant) for r and Tgi <^ for p. In the semiclassical 
limit we can therefore assume that r and p are statistically independent and the average 
consists of two independent integrals of r in the particle volume and of p on the Fermi surface 
\p \ = mvp- For a circle of radius R'm. d = 2 dimensions we obtain = \ and for a sphere 
of radius R'm. d = ?> dimensions the result is k = ^. One can also consider a (i-dimensional 

15 

ring [d = 1, 2, 3) on the xj/-plane with a thickness much smaller than R but larger than the 
mean free path. The integrations then yield = \- 

In the ballistic case we need the two constants k and k which appear only in the combi- 
nation K k. They are defined by 

^ ^ ('^\ m^R^yl = r dt {LM LM)q,v (CI) 
\Vp / Jo 

where the time dependence of the angular momentum is determined by the dynamics of 
the Hamiltonian without magnetic field. In Eq. ( |C1|) already the time reversal invariance 
was exploited. Let us now consider the classical ballistic motion of the electron in the 
quantum billiard. Let be the i-th reflection point on the surface (with z = 0, 1, 2, . . . and 
'^o = "f^it = 0)) and tj the corresponding time of the reflection. The classical path is then 
given by f{t) = fi + vptci where ti < t < tj+i and Cj is the unit vector in the direction of 
rj_i_i — Tj. We can now apply the semiclassical approximation to the average 



with 

L(ti) = mvp-. 

\ri+i - ri\ 

Let us now consider a surface that is very rough on a scale much smaller than the diameter 
of the cavity. In addition, we consider again a circle in two dimensions or a sphere in three 
dimensions. We assume that the surface-roughness implies that all reflection points are 
statistically independent and that the average with respect to fi is just an integration over 
a circle (or sphere) with radius R. Then in Eq. ( |(J2D only the term with i = yields a non 
vanishing contribution. Let d be the angle between ro and rl. In two dimensions, Eqs. (|C1|) 
and ( |(J2| ) yield after some obvious vector algebra 



j=0 



Kk = Uma§) cos=^(i?9)\ = — . (C3) 

\ ^2 ^ V2 Vaveragemd = 2 Syr ^ ' 

In three dimensions the angular momentum has also an x- and a y-component and we obtain 
Kk = ]- Um{U) cos'iU)) = ^ . (C4) 

3 \ ^ 2 ^ V 2 V average in d = 3 45 ^ ^ 

The relation (J4.2D has then the form 
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A-C.,^e) (C5) 



V Aoi? V$o 
with 

C2 = -J^ ^ 0.3676 , (C6) 
7r V 3 

C3 = W — ~ 0.2379 . (C7) 



The square root factor in Eq. (|C5|) can be calculated by standard expressions that relate 
the Fermi energy and the number n of electrons in the cavity with the result: 



I hvp 



n^l^ md = 2 , (C8) 
(6/7r)i/*^ md = 3 . (C9) 



The numerical results for an asymmetric stadium on page 176 of Ref. yield a relation 



A = 1.5 • ($/$o) for the first 500 energy levels. A comparison with ( |C5| ) and ( |C5| ) imply 
that the corresponding number of electrons is just n ~ 280 which is just in the middle of 
the considered energy range. The agreement is very good if we consider the simplicity of 
the above used "rough-surface" model. The numerical values ( P6| - |C7| ) should therefore be 
considered as (rather good) estimations for a typical quantum billiard. 
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FIGURE CAPTIONS 



Fig. 1 Conductance distribution p{T;t) for = 1 (Eq. ( ^lilD ) and t/tc = 0, 0.5, 1, 2, oo. 



Fig. 2a Transmission eigenvalue density Ri{T;t) for = 5 and t — oo (unitary limit) and 



the large A^-limit N/{'k^T{1 - T)) found in Ref. |U| . 
Fig. 2b The same as in Fig. 2a, but for t = (orthogonal limit). 

Fig. 3 Difference Ri{T] t) — Ri{T] oo) for = 5 and t/tc = 0, 0.5, 1, 2, oo. The curve with 
largest deviation corresponds to t/tc = 0. 

Fig. 4 2-point correlation function i?2(To,T;t) for Tq = 0.7 and 0.6 < T < 0.8 at = 5 
and t/tc = 10~^, 10~^, 10~^, 1, oo. One can see the beginning of the COE — > CUE 
crossover on small scales |T — To| when t/tc ^ 1- 

Fig. 5 The weal-localization suppression {g)cvE — {g(t)) ioi N = 1, 5, oo as a function of 

t/tc. 

Fig. 6 Conductance fluctuations {5g{tY) for A^ = 1, 5, oo as a function of t/tc- 
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